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ABSTRACT 

We present multi-frequency VLBA observations of two polarized Compact 
Symmetric Objects (CSOs), JOOOO+4054 and J1826+1831, and a polarized CSO 
candidate, J1915+6548. Using the wavelength-squared dependence of Fara- 
day rotation, we obtained rotation measures (RMs) of — 180±10 rad m~^ and 
1540±7 rad m~^ for the latter two sources. These are lower than what is ex- 
pected of CSOs (several 1000 rad m~^) and, depending on the path length of the 
Faraday screens, require magnetic fields from 0.03 to 6 fiG. These CSOs may be 
more heavily affected by Doppler boosting than their unpolarized counterparts, 
suggesting that a jet-axis orientation more inclined towards the line of sight is 
necessary to detect any polarization. This allows for low RMs if the polarized 
components are oriented away from the depolarizing circumnuclear torus. These 
observations also add a fourth epoch to the proper motion studies of JOOOO+4054 
and J1826+1831, constraining their kinematic age estimates to >610 yrs and 
2600±490 yrs, respectively. The morphology, spectrum, and component motions 
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of J1915+6548 are discussed in light of its new classification as a CSO candidate, 
and its angle to the line of sight (~50°) is determined from relativistic beaming 
arguments. 

Subject headings: galaxies: active - galaxies: ISM - galaxies: jets - galaxies: 
nuclei - radio continuum: galaxies 



1. Introduction 



Compact symmetric objects (CSOs) are a class of active galactic nuclei ( AGN) that have 



signif i cant visible jet or hots pot activity on either side of the central engine (j Conway et al. 



I994J : IWilkinson et al.l 11994) . These are typically < 1 kpc in size because they are young 
objects (<3000 v r: iPolatidis fc Conwav II2003I : iGugliucci et al.ll2005l ). In terms of the unified 
scheme of AGN (lAntonucci Ill993l ). there is little Doppler boosting of the jets and hotspots 
due to the orientation of the source with respect to the line of sight. This orientation also al- 
lows for studies of the hypothesized dust and gas torus that surrounds the central engine and 
lies perpendicular to the jet axis. The actual angle that the line of sight makes with the torus 
determines what features can be seen. Evidence for this torus comes from detections of broad 



HI absorption lines ([Taylor et al.lll999l : iPeck et al. 



2OO3I : I Gupta et al.l l2006l ) and free- free absorption ( 



2000 



Pihlstrom. Conwav. fc Vermeulen 



Peck. Taylor, fc Conway 111999^ towards 



CSOs. H2O megamasers also provide clues to circumnuclear tori in nearby AGN JLo 2005 
and references therein). 

This torus, and magnetic fields in the shocked, photoionized ISM around the lobes, 
may also give rise to the large Fa raday screens that depolarize emission from most CSOs 
(IBicknell. Dopita. fc O'Dea Ml997l ). The first detection of significa nt polarization in tw o 
CSOs in the COINS sample (CSOs Observed In the Northern Sky; IPeck fc TavloTlboOoh . 
JOOOO+4054 (2. 1%) and J1826-j-1831 (8.8%), occured on the side with the more prominent 
jet or hotspot (IGugliucci et al.l l2005l . hereafter GTPG). The opposite hotspot showed no 
polarization down to a 0.3 mJy limit. An y polarization that is detected sh ould have high 



rotation measures (several 1000s rad m~^; iBicknell. Dopita. fc O'Dea 1119971 . and references 
therein). In this paper, we attempt to calculate the Faraday rotation measures in these 
polarized components. 

Kinematic ages can be obtained by measuring the separation speed between hotspots 
over time, or the proper rnotion of a hotspot or jet component with respect to the core (e.g. 
Polatidis k Conwav II2003I : iNaeai et al.ll2006h . In GTPG, we confirmed ages for three CSOs 
between 130±47 and 3000±1490 yr, with another source appearing to have an age of 20±4 
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yr. The overall CSO age distribution seems to be disproportionately stacked towards the 
younger ages. Possible explanations are: there is a selection effect against older CSOs, the 
jet activity tends to d ie off after a certain period of time, o r CSOs have periods of deactiva- 



tion and reactivation ( Tingay. Edwards, fc Tzioumis 20031). Those that survive may evolve 



into Fanaroff-Riley II galaxies suc h as Cygnus A (IFanaroff fc Riley I Il974j : iReadhead et al. 



1996bl : iPolatidis fc Conway I l2003l ). In this paper, we refine the kinematic age estimates of 



JOOOO+4054 and J1826+1831 that were first calculated in GTPG. Also included are obser- 
vations of a new CSO candidate from the Second Caltech-Jodrell Bank Survey, J1915+6548, 
which also shows hotspot separation over three epochs and polarization in its more prominent 
hotspot. 

Throughout this discussion, we assume Ho=71 km s^^ Mpc~^, Qm = 0.27, and ^a= 
0.73. Linear sizes and velocities for sources with known redshifts have been calculated using 
E.L. Wright's cosmology calculator 0. 



2. Observations and Analysis 

Observations were centered on 4.8 GHz and 8.4 GHz on 12 February 2005 for JOOOO+4054 
and on 18 February 2005 for J1826+1831 using the VLBA§. Observations of J1915+6548 
were centered on 4.8 GHz, 8.4 GHz, 15.1 GHz, and 22.2 GHz with the VLBA on 11 Novem- 
ber 2004. Each frequency was separated into four IFs, and these IFs were paired such that 
the higher two frequencies were averaged during imaging as were the lower two frequencies, 
except at 22.2 GHz where all four IFs were averaged. Therefore, the frequencies used for the 
RMs in these observations were 4.6 GHz, 5.0 GHz, 8.2 GHz, 8.5 GHz, 14.9 GHz, 15.3 GHz, 
and 22.2 GHz. Observational parameters are presented in Table 1. 

Amplitude calibration of the data was derived from system temperatures and antenna 
gains. Fringe-fitting was performed with the AIPS task FRING on the calibrators OQ 208 
and 3C 84. D-term solutions were determined with the AIPS task LPCAL and the same 
calibrators OQ 208 and 3C 84. Absolute electric vector position angle (EVPA) calibra- 
tion was determined for JOOOO+4054 using the EVPAs of J1310+322 and for J1826+1831 
and J1915+6548 using the EVPAs of BL Lac listed in the VLA Monitoring Progra n@ 



^http:/ /www. astro. ucla.edu/'^wright/CosnioCalc. html 

^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation. 

http://www.vla.nr ao.edu/astro/calib/polar/ 
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flTavlor fc Mvers I l2000h . 
separately. 



Note that the EVPAs were corrected for each of the four IFs 



Results 



Faraday rotation was first noted by Michael Faraday when he passed polarize d light 
through a refractive medium in the presence of a magnetic field (IFaraday I Il933l ). The 
intrinsic polarization angle, xo, is observed as x such that 



'V. 



where A is the observed wavelength. The rotation measure, RM, is related to the electron 
density, ne, the net line of sight magnetic field in the environment, Sy, and the path length, 
dl, through the plasma, by the equation 



RM = 812 



ncB\\dl 



rad m ^ 



(2) 



where units are in cm~^, mG, an d pars ecs. A reasonable rig for radio galaxies is 10'^ cm~^ 
as estimated by I Zavala &: Taylor I (120031 ). Our lower limit for the path length is 0.3 pc, the 
approximate size of clumps within the Faraday screen in M87 (jZavala fc Taylor 1120021 ). Since 
the Faraday screer i for an AGN is now co nsidered to come from interactions of the jet with 
ambient material (IZavala fc Taylor II2004I ). a good upper limit for the path length is the jet 
radius, or 10 pc. Field strengths calculated with these parameters can be compared to the 
strength of a magnetic field in pressure balance with a thermal gas of the same and a 
temperature of 10^ K using 

— = UekT. (3) 

The 8.4 GHz images of JOOOO+4054 and J1826+1831 at full resolution are shown in 
Figures 1 and 2 with sticks representing polarization vectors. The 8.4 GHz images were then 
tapered to match the resolution of the 4.8 GHz images. The same circular beam was then 
applied to both images and spectral index maps were created. These are overlayed with 
the 4.8 GHz images in Figures 1 and 2. A plot for the rotation measure of the polarized 
component of J1826+1831 is shown in Figure 3. Source parameters for these CSOs and 
J1915+6548 are in Table 1. 



In GTPG, we attempted to obtain relative proper motions for CSOs in the COINS 
sample in order to calculate a kinematic age for each source. This assumes that the separation 
speed of the hotspots is uniform. This method is independent of the angle that the source 
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makes with respect to the hne of sight and of the distance to the source. This method also 
produces resuhs if the hotspot or outer jet component can be seen moving away from the core. 
Using 8.4 GHz data from three epochs spread across a five year period, we calculated the 
kinematic ages of three CSOs with reasonable certainty, while the rest provided lower limits. 
These new observations provided us with a fourth epoch for JOOOO+4054 and J1826+1831 
at 8.4 GHz, extending the time baseline to seven years. This provides better estimates and 
limits of the ages of these two CSOs. 

Models for JOOOO+4054 and J1826+1831 were fit to the visibility data for the March 
2000 epoch, since this was close to the middle of the time baseline and had a high dynamic 
range. These models were comprised of elliptical Gaussians. Visibility data from each epoch 
were fit to a model with the same size and shape ellipses so that only the positions and fluxes 
of these were allowed to vary with time. Uncertainties in position for each component were 
calculated from the signal-to-noise ratios and the synthesized beam. The positions of these 
components with respect to a reference component were fit with a weighted least squares line 
where the slope of the line provides the relative speed. Motions are considered significant 
if they are at least Sa above the errors and if the reduced chi-squared is nearly 1. Errors 
are dependent on the image noise and the individual flux of each component. If large errors 
are present, then the motions are considered upper limits, so that they give rise to ages that 
are lower limits. Table 3 includes the modelfit parameters, and Figure 4 gives plots of the 
proper motion of two components in J1826+1831. 

Total intensity images of J1915+6548, the new CSO candidate, are presented in Figure 5 
with sticks representing polarization vectors overlayed. The integrated total power spectrum, 
as well as that of component A alone, is shown in Figure 6. Since the brightest hotspot was 
polarized at 8.4 GHz, 15.1 GHz, and 22.2 GHz, a rotation measure plot is given in Figure 
7. Figure 8 presents the proper motions of hotspot separation in this source at 4.8 GHz 
over an 11 year timespan. Table 4 gives hotspot brightnesses for the polarized sources and 
the unpolarized sources of the COINS sample for comparison in §4. Plots of /3 vs. 6 for the 
separation speeds and orientations of J1826+1831 and J1915+6548 are presented in Figure 9. 



3.1. JOOOO+4054 



This CSO was identified as such in lDallacasa et al.l (120021 ) and GTPG, and polarization 



was detected at 8.4 GHz. It is associated with a galaxy of magnitude 21.4 (jStickel fc Kuehr 



19961 ). The core has still not been positively identified from the spectral index map (see Fig. 
1). From their compact morphologies, either Bl or B2 could be the core, but their spectral 
indices are rather steep (a ~ —0.6; Si, oc z^"). However, there appears to be a flattening 
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a ~ —0.4) between components B l and B2, so the core may be located there (see inset 



Fig. l.b). In lDallacasa et aLl (|2002| ). the combined flux of Bl and B2 (component Ce in that 
paper) is 71 mJy at 1.6 GHz. When compared with the fluxes of Bl and B2 here at 8.4 GHz, 
this yields a spectral index for that region of a ~ —0.25. High dynamic range imaging at 
15 GHz should provide a positive identiflcation. 

Polarization was detected in the southern hotspot at 8.4 GHz with an intensity of 
2.1 mJy (see Fig. 1). This is 1.2% of the intensity of the hotspot. There is no detectable 
polarization at 4.8 GHz down to the 3(j level (~ 0.2 mJy, or 0.06%). The 4.8 GHz image also 
shows a strange morphology in the southern lobe. There appears to be a hole, or depression, 
in the emission below the bright hotspot. If this is not an a rtifac t in the data, it resembles the 



region in the eastern lobe of 4C 31.04 by lGiroletti et al.l (120031 ) . where it is speculated that 



such a hole could be created by a dense molecular gas that is impenetrable to radio emission 
or to the plasma itself. However, the morphology may also be the effect of edge-brightening 
at component C if the plasma is running into denser material. This material may also be 
deflecting the jet to the east, as seen in the 8.4 GHz image. Polarization was only detected 
at 8.4 GHz, so the RM could not be determined and the true orientation of the magnetic 
field is not known. If the RM was known to be as low as in J1826+1831 (—180 rad m~^, see 
§3.2), for example, then there would be little rotation from the original polarization angle 
to 8.4 GHz, so the magnetic field would lie roughly perpendicular to the electric vectors in 
Figure la. Then, the detectable polarization in this region would be due to a compression 
of the ambient magnetic field at a shock front where the jet is colliding with the denser gas. 
Future observations at 15 GHz should allow us to determine the rotation measure for this 
region. 

The fourth epoch of 8.4 GHz data in the proper motion study of this object did not 
provide a good fit for the separation speed of the hotspots. However, there is an upper 
limit of 0.066 mas yr~^ of component A away from component C. This speed limit and a 
hotspot distance of 40.33 mas provides a lower limit for the kinematic age of 610 yrs. This 
is higher than was previously estimated (280 yrs in GTPG). Although the redshift is not 
known, a reasonable estimate of z ~ 0.5 can be made for the typical redshift of a CSO host 



(lAugusto et al.ll2006l ). If this is the case, the projected distance between components A and 



C is 242 pc and A is moving with speed v < 2.0c. 



3.2. J1826+1831 



This is the most significantly polarized of these three objects. Component C, most likely 
a jet component, has a polarized intensity of 2.3 mJy at 8.4 GHz, or 8.5% of its fiux (see 
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Fig. 2). This component has a polarized intensity of 1.4 mJy at 4.8 GHz, making it 3.7% 
polarized. A least-squares fit for the rotation measure of — 180±10 rad m~^ at the peak of 
component C is p resented in Figure 3. T his is not an uncommon rotation measure for a 



typical quasar jet (IZavala fc Taylor II2004I ). The pair of angles at 8.2 and 8.5 GHz suggest 
a higher rotation measure, but this requires for a number of 180 degree turns to be put in 
between 4.6 and 5.0 GHz and between 5.0 GHz and 8.2 GHz. These turns can be introduced 
because the polarization vectors give an orientation, but not a direction, of the electric field 
such that X ^iid x±180 are indistinguishable. Depending on the number of turns used in 
this data, rotation measures as high as 6000 rad m~^ are plausible. One must take caution 
with this, however, since any number of turns can be introduced to provide a false good 
fit. Therefore, — 180±10 rad m~^ is a conservative estimate. Measurement of polarization 
at another frequency may help in determining the correct rotation measure. It may also be 
true, however, that we only see polarized emission from this component because it is free of 
the high rotation measures that may be present closer to the core. 

If a path length of 0.3 pc is assumed with Ue = 10^ cm~^, the magnetic field strength is 
0.7 /iG, using Eqn. 2. With the upper limit of the path length, 10 pc, the field strength is 
only 0.03 //G. Both of these are much smaller than the strength of a magnetic field that is 
in pressure balance with a 10^ K gas of the same electron density, ~200 /iG, using Eqn. 3. 

A spectral index map was made between 4.8 and 8.4 GHz which is presented in Figure 
2. This shows the core, B, to be a flat spectrum component (a ~ 0.6), while the rest of the 
jet and counterjet are steep spectrum (a ^ —0.6 to —1.7). 

The proper motions of components C and D with respect to B were measured using the 
four epochs of 8.4 GHz data available to us. D was found to be moving at 0.015±0.003 mas yr~^. 
With a projected separation between B and D of 41.87 mas, this corresponds to a kinematic 
age for the CSO of 2600±490 yr. This refines our earlier estimate of 3000±1490 yr. The jet 
component C is moving along at a faster rate of 0.032±0.003 mas yr~^ (see Fig. 4). The 
modelfitting of components works best for bright, compact components, of which component 
A is neither, so no proper motions can be fit. If we assume again that z ^ 0.5, D is moving 
away from B at ~0.45c and is 255 pc away. C is moving away from B at ~0.95c. However, 
this source has no detected optical counterpart, so its redshift could be much higher. 



3.3. J1915+6548 



This new CSO candi date was first listed in a survey with the 300-ft Green Bank 
Telescope at 4.8 GHz by iBecker. White, fc Edwards I (Il99ll ) It was observed in the Sec- 
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ond Caltech-Jod r ell B ank Survey (CJ2) which presented a spectrum of total intensity in 



Henstock et al. I (I19951). Its host is a Se yfert 1 galaxy with a magnitude of 18.2 and red- 



shift of z = 0.486 (IHenstock et al. 1119971 ). New observations were taken with the VLB A to 



investigate the non-linear morphology in the CJ2 image. 

We identify Component A as a steep spectrum hotspot (a ~ —0.9; See Fig. 6) that 
dominates the flux of this object at all four frequencies. Components B and C are steep 
spectrum {a ~ —1.2 and —0.9, respectively), and we propose that they are also part of the 
eastern lobe. Component D is very steep spectrum (a ^ —1-6), and, based on morphology, 
we propose it to be the counter-hotspot to A. There is no evidence of a compact, fiat 
or inverted spectrum core, but there are a growing number of CSOs that have two steep 
spectrum hotspots but no visible core (e.g. J0620-f2102, JlllH-1955; GTPG). There is also 
a trail of extended emission between A and D at 4.8 and 8.4 GHz that hints at a jet-like 
structure, as well as emission further east of A in the 4.8 GHz image which we interpret as 
a sign of earlier activity (see Fig. 5). Although the morphology is similar to that of a CSO, 
the hotspot intensity ratio at 4.8 GHz is 12:1 and is 26:1 at 8.4 GHz. T hese fall short of the 
10:1 criteria set for CSOs in the COINS sample jPeck fc Tavlor Iboooh . 



To determine the synchrotron age of the source, we present the total intensity spectrum 
of J1915+6548 in Figure 6 using total intensities from this paper and data points from 
the NASA/IPAC Extragalactic Database (NED)0. Using the break fr equency, 17 G Hz, the 
minimum energy magnetic field of the plasma can be calculated as in iMiley I (Il980l ) using 



Bri 



1.4 X 10-^(1 + z 



1.1, ,0.22 







OxOyS 



2/7 



(4) 



where B is in gauss, uq is in GHz, Fq is in Jy, 6^ and 6y are the dimensions of an elliptical 
component in arcseconds, and s is t he path length in kp c. With this, the synchrotron age 
of the source can be calculated as in iMurgia et al.l (119991 ) with 

1610S°-5 



[i?2 + S^,,^]K(l + .)]V2 (5) 

with B in /iC, iyf,r in GHz, and Bqmb = 3.25(1+2;)^. At 4.8 GHz, the eastern hotspot is 
0.291 Jy and can be approximated by an ellipse 0.01 mas by 0.006 mas. We assume the 
path length to be equivalent to 0.01 mas, or 0.07 kpc at a redshift of 0.486. This yields a 
magnetic field of ~7 mC and a synchrotron age of ~540 yrs. 

Component A is polarized at 8.4 GHz (3.0 mJy or 1.6%), 15.1 GHz (2.3 mJy or 2.8%), 
and 22.2 GHz (1.8 mJy or 4.0%). The polarized fiux is too low to be reliable at 4.8 GHz 



''NED is operated by the Jet Propulsion Laboratory, Caltech, under contract with NASA. 
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(1.3 mJy or 0.5%). A plot of the rotation measure is given in Figure 7. The rotation 
measure of 1540±7 rad m~^ comes from a direct measurement of the EVPAs at the center of 
the polarized component without any introduced 180 degree wraps. This is more consistent 
with the high RMs expected from a CSO. With a path length of 0.3 pc and the same electron 
density as above, this corresponds to a magnetic field strength of 6 fxG. At the upper limit 
for path length of 10 pc, the magnetic field strength is 0.2 /iG. These are still much lower 
than what is expected in the central regions of a radio galaxy (200 /iG, see §3.2). 

The visibility data at 4.8 GHz was modelfit at this epoch and at tw' o othe r epochs. 



1993.444 and 1995.691, the first of which was presented in iHenstock et al. I (119951 ). Proper 
motions were fit with a least squares line with a slope of 0.035±0.004 mas yr~^. With 
a redshift of z =0.486, this apparent speed corresponds to 1.02±0.11 c. Although this is 
higher than typical hotspot separation speeds (GTPG), it yields a good fit (see Fig. 8) and 
a kinematic age of 940±110 yr. 



4. Discussion 
4.1. CSO Orientations 



Depolarization in CSOs is consistent with the lBicknell. Dopita. fc O'Deal (119971 ) model. 
Variations in rotation measure across the lobe can cause this depolarization, and this may 
be due to magnetic field reversals in the post-shock ISM. The obscuring torus may also play 
a role in depolarization, and this would be dependent on torus scale height, opening angle, 
and orientation. We rule out depolarization by a thin disk of material, with a scale height 
much less than the radio source size, since we frequently see that both sides are depolarized, 
whereas a thin disk would only cover one side. If a Faraday screen is to produce the observed 
polarization asymmetry in these sources, then it requires a scale height of about half the total 
source size. More detections of polarization in CSOs are needed to constrain its geometry. 

The existence of polarized emission in these three sources sets them apart from the 
rest of the COINS sample. However, this is not the only difference. In GTPG, we noted 
that JOOOO+4054 and J1826+1831 appeared to have greater hotspot intensity ratios than 
their unpolarized counterparts. Since the polarization is detected on the side of the more 
prominent hotspot, this can be attributed to a shorter path length through the circumnuclear 
torus and, consequently, a lower Faraday depth. Table 4 lists the brightnesses of each hotspot 
and the ratio of the dimmer one to the brighter one for unpolarized and polarized CSOs. 
There is a bimodal distribution in that the polarized sources presented in this paper have 
ratios <0.1 and the unpolarized COINS have ratios >0.3. Although the sample size for 
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polarized CSOs is small, this is consistent with the observed polarization being due to an 
orientation effect. As the angle between the jet axis and line of sight decreases, the forward 
hotspot may be more Doppler boosted and the receding hotspot Doppler dimmed, assuming 
mildly relativistic bulk motions. So, the leading hotspots in the three polarized sources are 
oriented such that they are relatively free from Faraday depolarization by the torus that 
surrounds the central region. However, to explain the depolarization over many viewing 
angles in the other CSOs, a large torus is needed to extend over hundreds of parsecs of 
jets and l obe. HI absorption has been s een at this distance from the central engine (e.g. 



1946+708 IPeck. Tavlor. fc Conwav \\l99% . 



In order to quantify depolarization in CSOs, we look at Faraday beam depolariza- 
tion. We can calculate the necessary rotation measure gradient across the synthesized beam 
in order to get a rotation of one radian. For the unpolarized CSOs in the COINS sam- 
ple, with an average beamwidth of 1.5 mas at 8.4 GHz, the RM gradient need only be 
~500 rad m~^ mas~^. So, the RM of J1915+6548 may be reasonable for a CSO, whereas the 
low rotation measure of J1826+1831 is still puzzling. An RM of 500 rad m~^ corresponds 
to a magnetic field strength of 2 /iC with a path length of 0.3 pc and 0.06 ^G with 10 pc. 
These are still much lower than the strength of a magnetic field in pressure balance with a 
10^ K gas. 

We use the assumption that relativistic beaming is in effect to quantify the orientations 
of these sources. We can match up what appear to be jet components from the approaching 
and receding sides of the core and compare their fluxes {Sa and Sr) as 

Sa /l + /3cos6' 



C \, n n, (6) 
br Vl ~ PCOS 

where (3 is the space velocity, 9 is the angle to the line of sight, a is the spectral index, 
and n is either 2 or 3. Models for continuous jets are best fit by n = 2 and jets of discrete 
components by n = 3. One can use 

fla _ da _ 1 + COS 6' 
Hr dr 1 — COS 9 ' 

where /i is the apparent motion and d is the distance from the core. However, this requires 
a well known position for the center of activity. One can also use the hotspot separation, 

A^sep = \l^a\ + ll^rli SUch that 

/i \ 2/9 sin ^ 

Vsep = l^sepDail + z) = — , 8 

1 — cos^ 9 

where Da is the angular size di stance to the source, z is the redshift, and Vsep is the angular 



separation speed in units of c ( iTaylor fc Vermeulen 1119971 ). These can provide constraints 
on (3 and 9. 
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We know neither the core location nor the redshift for JOOOO+4054. So, only Eqn. 6 can 
be applied, using the fluxes of the hotspots as determined by the model in Table 3. A locus 
of P and 6 is plotted in Figure 9a. Note that a lower limit is plotted when n = 2 and an 
upper limit when n = 3. The errors in the measured fluxes are negligible compared to this 
spread. This shows that, approximately, /3 > 0.3 for all values of 6. 

For J1826+1831, we used the brightnesses of each hotspot, components A and D, with 
Eqn. 6 as well as their distances from the core with Eqn. 7 to calculate loci of /3 and 6 
(Figure 9b). Note that these two equations have a similar functional form, so no tight 
constraints can be determined. This gives an upper limit for the angle to the line of sight 
of ~75°. The true space velocity wil l be >0.3 c for any angle. T his is higher than typical 



values found for othe r CSOs (^O.lc lOwsianik &: Conway I Il998l ) and other radio galaxies 



(<0.1c IScheuer I Il995l ). If the value for fi^^p of component D (0.015 mas yr~^)is used in 
Eqn. 8 with a typical CSO redshift {z ^ 0.5), the projected separation velocity is ~0.4c. 
This is consistent with our results. J1826+1831 cannot be very close (for example, z = 0.01), 
for then the separation velocity would be unusually low and inconsistent with Figure 9b. 

We used the apparent hotspot separation speed of 1.02±0.11c for J1915+6548, z =0.486, 
and Eqn. 8 to plot the locus for /3 and 9 in Figure 9c along with the loci for the two brightness 
models with Eqn. 6. These two sets of curves intersect when /3 ~ 0.55 and 9 ^ 50°. Such a 
high space velocity is uncharacteristic of lobes. However, it is possible that the hotspot being 
measured is a transient feature, and therefore moving more quickly, whereas the lobe as a 
whole is moving at a much slower speed. This can occur if the jet that is fee ding the hotspo t 



has changed its orientation slightly and is drilling out a new part of the lobe (IScheuer Ill995l ). 
Alternatively, we may be measuring the separation of jet components rather than hotspots 
as there may be no visible hotspots. This example shows how the orientation angle of a CSO 
can be measured if mild relativistic beaming is a factor. 



4.2. CSO Environments 



An alternate th eory for the sma ll sizes of CSOs is that they are old sources frustrated 
by a dense medium (jCarvalho Ill994j ). This dense m edium may in fa ct be asymmetric, pos- 
sibly as a result of a galaxy merger or interaction (jCarvalho 1 119981 ). In an extreme case, 
the ages of these sources would be comparable to the ages of larger, classic radio doubles. 
Readhead et al. I (Il996al ) pointed out that if this were t he case, CSOs would have more 
spherical morphologies. Also, hotspot separation speeds in iPolatidis &: Conway I (120031 ) and 
GTPG show that the hotspots are still moving too quickly to be confined by a dense medium 
and be older than a few thousand years. However, this does not rule out high densities or 
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asymm e tries in the medium around CSOs as evidence by lOrienti. Dallacasa. fc Stanghellini 



( 120061 ) ■ iBicknell. Dopita. fc O'Dea I (119971 ) propose that CSOs are frustrated but not con- 
fined by the interstellar medium (ISM). Their model assumes a dentist-drill explanation for 
jet-lobe interactions and predicts low polarization. This is because the ionized gas surround- 
ing the lobes in this model produces large variations in Faraday rotation measures across 
the source, thus depolarizing the emission. Interstellar magnetic fields can play an integral 
role if there are a large number of magnetic field reversals across the source which would 
produce a varied RM structure in our maps. Since the polarization is spatially isolated in 
these CSOs, no such maps can be made. 

Asymmetries in the CSO environment may cause one hotspot to be considerably more 
polarized in these few sources by interactions with a dense ISM. If the EVPAs for J1915+6548 
are extrapolated back to zero wavelength, the observed angle of the electric field is —88°. 
So, the orientation of the magnetic field in the image would be nearly north-south, which 
is approximately perpendicular to the source orientation. This may be an indication of a 
collision between the bright hotspot and a dense medium that orders the magnetic field 
along the axis of compression. This is also observed in J1826+1831 where the EVPA at zero 
wavelength is 75°, so the magnetic field is oriented at —25°. This is also nearly perpendicular 
to the source orientation, but this effect is for the bright jet component, not the hotspot. 
Components may also be brighter (by Doppler boosting) if they are moving more quickly 
through a less dense ISM. However, the indistinguishable motions of the polarized component 
in JOOOO+4054 seems to refute that. 

This asymmetry in environments does not, however, explain the bimodal distribution 
in Table 4, that is, the fact that the ratios of the dimmer hotspots to the brighter hotspots 
of each source are much lower for polarized CSOs than for unpolarized CSOs. This is more 
naturally explained by Doppler boosting. Also, the relativistic hotspot speeds of J1915-I-6548 
(see Fig. 9) could not exist in a dense enviro nment. Finally, we do not see rotation measures 



in the several thousands as predicted by the iBicknell. Dopita. fc O'Dea I (119971 ) model. 



5. Conclusions 

In this study we have investigated the polarization properties of two CSOs, JOOOO-l-4054 
and J1826-I-1831, and a CSO candidate J1915-I-6548. The Faraday rotation measures that 
were observed for J1826+1831 and J1915-I-6548 are lower than what was previously expected 
for CSOs, — 180±10 rad m~^ and 1540±7 rad m~^, compared to predicted values of several 
thousand rad m~^. These imply low external magnetic field strengths, 0.03 to 6 /iC, de- 
pending on what estimates for electron density and path length are used. However, only 
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about 500 rad m ^ mas ^ is needed for beam depolarization to be a major factor for most 
unpolarized CSOs. 

Since these polarized sources have significantly more asymmetric hotspot brightnesses 
than their unpolarized counterparts, with ratios of <0.1 as opposed to >0.3, it is plausible 
that Doppler boosting is in effect for the brighter hotspots. Then, their jet axis orientations 
arc closer to the line of sight than for other CSOs such that the polarized hotspot is free of 
obscuration from the circumnuclear torus. We put contraints on (3 and 9 for all three sources 
(Figure 9) such that for JOOOO+4054 and J1826+1831, P > 0.3. More information, specif- 
ically the redshift, must be known in order to make tighter constraints. For J1915-I-6548, 
P fa 0.5 and ^ 50°. 

We extended the time baseline for proper motions of the two objects from the COINS 
sample from five to seven years. This provided a better age estimate for J1826+1831 of 
2600±490 yrs and a better lower limit for the age of JOOOO-l-4054 at 610 yrs. Extending 
the time basehne for the proper motions of a larger sample of CSOs will reduce some of the 
uncertainties and lower limits of their ages so that the distribution of ages in the sample can 
be better determined. 

J191 5+6548 was determined to be a CSO candidate because of the presence of symmet- 
ric, steep spectrum hotspots and the lack of a compact, fiat or inverted spectrum core. If 
confirmed, this will be the third CSO with significant polarization. Although spectral age 
arguments and proper motions from an eleven year time baseline suggest an age between 
approximately 600 and 1000 yrs, the hotspot separation speed is most certainly relativistic 
and Doppler boosting is occuring in the brighter hotspot. Our age estimates may be under- 
estimating the true source age if the brightest hotspot is undergoing renewed activity from 
a slowly moving jet or if it is encountering a particularly dense medium, reaccelerating the 
electrons and giving a lower synchrotron age. If it were the case that the hotspots are tran- 
sient features and therefore appear younger and faster moving than the radio lobe, then the 
kinematic ages of many CSOs would be underestimated. However, with such small sources, 
it is possible that a true classical radio lobe has not yet been created and the hotspots are 
the only indicator of the current phase of activity in CSOs. 

NEC gratefully acknowledges support from the NRAO Graduate Summer Student Re- 
search Assistantship. ABP thanks NRAO for travel support and hospitality during part of 
this project. This research has made use of the NASA Astrophysics Data System. 
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Fig. 1. — Total intensity contours of JOOOO+4054. (a) Electric polarization vectors over 
8.4 GHz contours. A vector length of 1 mas corresponds to a polarized flux density of 0.05 
mJy beam~^; the smallest ticks correspond to ~0.6 mJy beam~^. Contour levels begin 
at 0.25 mJy beam~^ and increase by factors of 2. (b) Spectral index map over 4.8 GHz 
contours where oc u". Contour levels begin at 0.7 mJy beam~^ and increase by factors of 
2. The inset rescales the spectral indices for the Bl and B2 components. A star indicates 
the reference component for motions. 
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Fig. 2. — Total intensity contours of J1826+1831. (a) Electric polarization vectors over 
8.4 GHz contours. A vector length of 1 mas corresponds to a polarized flux density of 0.10 
mJy beam~^; the smallest ticks correspond to ~0.8 mJy beam~^. Contour levels begin at 
0.25 mJy beam~^ and increase by factors of 2. (b) Spectral index map over 4.8 GHz contours 
where oc z/". Contour levels begin at 0.6 mJy beam~^ and increase by factors of 2. A star 
indicates the reference component for motions. 
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Fig. 3. — Plot of polarization angle versus wavelength squared for component C of 
J1826+1831 where the slope of the least squares line is the rotation measure. 





Fig. 4. — Proper motion plots of J1826+1831 with four epochs of 8.4 GHz models to the 
visibility data, (a) The slope for component D along the x axis is — 0.014±0.002 mas yr~^ 
with a reduced chi squared of 0.912. The slope along the y axis is 0.008±0.004 mas yr~^ 
with a reduced chi squared of 7.487. (b) The slope for component C along the x axis is 
— 0.031±0.003 mas yr~^ with a reduced chi squared of 2.385. The slope along the y axis is 
0.005±0.008 mas yr~^ with a reduced chi squared of 17.193. 





Fig. 5.— Total intensity plots of J1915+6548 at (a) 4.8 GHz, (b) 8.4 GHz, (c) 15.1 GHz, 
and (d) 22.2 GHz with electric polarization vectors overlayed. Contour levels begin at 
0.4 mJy beam~^ in (a) and (b) and at 5.5 mJy beam~^ in (c) and (d) and increase by factors 
of 2. A vector length of 1 mas corresponds to a polarized flux density of (b) 0.42 mJy bcam"-'^, 
(c) 0.83 mJy beam~^, and (d) 1.7 mJy beam~^. The bar in the lower right represents 5 mas 
or ~45 pc. 
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Fig. 6. — Total power spectrum of J1915+6548 using data points from the NASA Extra- 
galactic Database and the data presented in this paper (filled squares). The break frequency 
is 17 GHz. Overlayed is the spectrum of component A which has a spectral index of —0.9 
(open squares). 
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Fig. 7. — Plot of polarization angle versus wavelength squared for component A of 
J1915+6548 where the slope of the least squares line is the rotation measure. Note that 
X and y scales are different. 
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Fig. 8. — Proper motion plots of J1915+6548 with three epochs of model fit to the 4.8 GHz 
visibility data. The slope for component D along the x axis is — 0.024±0.002 mas yr^^ with 
a reduced chi squared of 0.118 The slope along the y axis is — 0.026±0.005 mas yr~^ with a 
reduced chi squared of 1.153. 



-25- 
JOOOO+4054 




20 



40 



60 



80 



100 



en 

J1826+1831 



1.0 



0.8 



0.6 - 



u 
> 

^ 0.4 



0.2 - 



0.0 




(b) 



Arm Length Ratio 
Brightness Ratio 



20 



40 



60 



80 



100 



en 

J1915+6548 



1.0 



0.8 



0.6 - 



u 
> 

"■^^ 0.4 



0.2 - 



0.0 




22^ 



(c) 



Motion Ratio 
Brightness Ratio 



20 



40 



60 



80 



100 



Fig. 9. — Plots of loci for intrinsic jet velocity, (3, vs. angle to the line of sight, 6, for 
(a) JOOOO+4054, (b) J1826+1831, and (c) J1915+46548. The brightness ratios are calculated 
from two different models with the hatching to indicate that the values lie between the lines. 
The arm length and motion ratios are hatched in between the upper and lower limits due to 
experimental error. Cross-hatching indicates where the regions overlap. Single-hatching is 
removed from (b) for clarity. 
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Table 1. Observational Parameters 



Source 


Date 


Preq. 


Time 


Bandwidth 


rms 






(GHz) 


(min) 


(MHz) 


(mJy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


JOOOO+4054 


20050212 


4.6 


200 


16 


0.18 




20050212 


5.0 


200 


16 


0.12 




20050212 


8.2 


200 


16 


0.12 




20050212 


8.5 


200 


16 


0.09 


J1826+1831 


20050218 


4.6 


212 


16 


0.17 




20050218 


5.0 


212 


16 


0.18 




20050218 


8.2 


212 


16 


0.20 




20050218 


8.5 


212 


16 


0.19 


J1915+6548 


20041111 


4.6 


58 


16 


0.23 




20041111 


5.0 


58 


16 


0.18 




20041111 


8.2 


39 


16 


0.17 




20041111 


8.5 


39 


16 


0.15 




20041111 


14.9 


97 


16 


0.26 




20041111 


15.3 


97 


16 


0.21 




20041111 


22.2 


135 


32 


0.18 



*Notes - (1) J2000 source name; (2) Date of observation; (3) 
Frequency in GHz; (4) Integration time in minutes; (5) Bandwidth 
in MHz; (6) rms noise in mJy. Note that the rms for each fully 
averaged frequency, 4.8, 8.4, and 15.1 GHz, is generally a factor of 
\f2 lower. 
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Table 2. Source Parameters 



Name 


RA 


Dec 


ID 




z 


'^5GHz 


%GHz 


'^15GHz 


'^22GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


JOOOO+4054 


00 00 53.081551 


+40 54 01.79335 


G 


21.4 




521 


322 






J1826+1831 


18 26 17.710882 


+ 18 31 52.88973 








427 


279 






J1915+6548 


19 15 23.819114 


+65 48 46.38505 


G 


18.2 


0.486 


331 


202 


104 


54 



*Notes - (1) J2000 source name; (2) Right ascension and (3) Declination in J2000 coordinates from the VLBA 
Calibrator Survey by Beasley et al. 2002; (4) Optical host galaxy identification; (5) Optical magnitude; (6) Redshift; 
(7) Total flux density at 4.8 GHz in mJy; (8) Total flux density at 8.4 GHz in mJy; (9) Total flux density at 15.1 GHz 
in mJy; (10) Total flux density at 22.2 GHz in mJy. 
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Table 3. CSO Model Parameters 













S 


P 




V 


Kinetic Age 


Source 


Component 


(mas) 


(mas) 


(deg.) 


(mJy) 


(mJy) 


(mas yr~^) 


(c) 


(yr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


JOOOO+4054 


A 


2.28 


1.71 


-3.4 


14 


<0.2 


<0.066 




>610 




Bl 


0.91 


0.42 


-83.6 


19 


<0.2 










B2 


0.63 


0.63 




28 


<0.2 










C 


3.01 


2.41 


-65.4 


205 


2.1 


Reference 






J1826+1831 


A 


8.05 


4.51 


24.6 


12 


<0.5 










B 


0.45 


0.45 




11" 


<0.5 


Reference 








C 


0.96 


0.77 


43.3 


32 


2.3 


0.032±0.003 




450±43'' 




D 


1.74 


1.11 


-71.3 


115 


<0.5 


0.015±0.003 




2600±490 


J1915+6548 


A 


0.46 


0.32 


25.2 


183 


<0.3 


Reference 








B 


1.88 


0.86 


29.8 


72 


<0.3 










C 


1.65 


1.37 


24.9 


36 


<0.3 










D 


1.30 


0.88 


17.6 


15 


<0.3 


0.035±0.004 


1.02±0.11 


940±110 



*Notes - (1) J2000 source name; (2) Component name; (3) Major and (4) minor axes of Gaussian model com- 
ponent; (5) Position angle of major axis; (6) Integrated flux density of Gaussian model component; (7) Polarized 
intensity, or 3a limit; (8) Relative proper motion: (9) Relative proper motion in terms of c if z is available; (10) 
Kinematic age estimate. For JOOOO+4054 and J1826+1831, fluxes are at 8.4 GHz at 2005.118. For J1915+6548, 
fluxes are at 4.8 GHz at 2004.863. 

"^Correction: Gugliucci et al. (2005) erroneously listed the flux densities of component B as 80, 90 and 70 mJy in 
epochs 1, 2, and 3. The correct fluxes are 8, 9, and 7 mJy, respectively. This gives a core fraction of 3%. 

^This is not a source age, but a kinematic age estimate for a jet component. 
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Table 4. Hotspot Ratios for Polarized and Unpolarized CSOs 







c 


c 


-tiaiio 


CyOrc 




Source 


(mJy) 


(mJy) 


(S2\ 


Fraction 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Polarized Sources 


JOOOO+4054 


205 


14 


0.07 


<0.0006 




J1826+1831 


115 


12 


0.10 


0.04 




J1915+6548 


153 


6 


0.04 


<0.001 


Unpolarized Sources 


J0003+4807 


38 


13 


0.34 


0.04 




J0204+0903 


120 


68 


0.57 


0.20 




J0427+4133 


65 


28 


0.43 


0.86 




J0620+2102 


156 


88 


0.56 


<0.001 




J0754+5324 


39 


36 


0.92 


<0.003 




Jllll+1955 


98 


76 


0.78 


<0.002 




J1143+1834 


130 


99 


0.76 


<0.001 




J1414+4554 


64 


51 


0.80 


<0.003 




J1546+0026 


213 


110 


0.52 


0.40 




J1734+0926 


183 


119 


0.65 


<0.0006 




J1816+3457 


135 


59 


0.44 


<0.001 




J2203+1007 


145 


53 


0.37 


<0.001 



*Notes - (1) Category; (2) J2000 source name; (3) Integrated flux density 
of Gaussian model component of brighter hotspot at 8.4 GHz; (4) Integrated 
flux density of Gaussian model component of dimmer hotspot at 8.4 GHz; 
(5) Ratio of dimmer hotspot as compared to brighter hotspot; (6) Fraction 
of total flux that is attributed to the core (using 3a limit for core flux if not 
detected). Fluxes for unpolarized CSOs from 2000.227 in GTPG. Fluxes 
for polarized CSOs and candidate from this paper. 



